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EVALUATION 


This  research  has  quantitatively  established  the  validity  of  adaptive 
excision  of  additive  narrowband  interference  from  direct  sequence  spread 
spectrum  communication  signals.  More  importantly,  it  has  shown  that  this 
strategy  is  suitable  for  dispersive  channels  such  as  the  skywave  High 
Frequency  (HF)  radio  communication  channel.  Consequently,  this  work  is 
directly  applicable  to  the  accomplishment  of  RADC  TPO  in  the  area  of 
sub-UHF  communications,  especially  HF  radio  communications.  We  anticipate 
that  the  research  described  in  this  report  will  be  directly  applied  in 
our  forthcoming  exploratory  development  work  in  HF  signal  processing  and 
in  our  advanced  development  new  initiative  in  HF  communications. 


/(JOHN  T.  GAMBLE 
Project  Engineer 
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I.  INTRODUCI  rON 

Spread  spectrum,  direct  sequence  or  pseudo-noise  (PN)  modulation  is 
employed  in  digital  communication  systems  to  reduce  the  effects  of  inter¬ 
ference  due  to  other  users  and  intentional  jamming.  When  the  interference 
is  narrow  hand  the  cross-correlation  of  the  received  signal  with  the 
replica  of  the  PN  code  sequence  reduces  the  level  of  the  interference  by 
spreading  it  across  the  frequency  band  occupied  by  the  PN  signal.  Thus, 
the  interference  is  rendered  equivalent  to  a  lower  level  noise  with  a 
relatively  flat  spectrum.  Simultaneously,  the  cross-correlation  operation 
collapses  the  desired  signal  to  the  bandwidth  occupied  by  the  information 
signal  prior  to  spreading. 

The  interference  immunity  of  a  PN  spread  spectrum  communication  sys¬ 
tem  corrupted  by  narrow  band  interference  can  be  further  improved  by 
filtering  the  signal  prior  to  cross -correlation,  where  the  objective  is 
to  reduce  the  level  of  the  interference  at  the  expense  of  introducing 
some  distortion  on  the  desired  signal.  This  filtering  can  be  accomplished 
by  exploiting  the  wideband  spectral  characteristics  of  the  desired  PN 
signal  and  the  narrow  band  characteristic  of  the  inteference.  Since  the 
spectrum  of  the  PN  signal  is  relatively  flat  across  the  signal  frequency 
band,  the  presence  of  a  strong  narrow  band  interference  is  easily  recog¬ 
nized.  Then,  the  interference  can  be  suppressed  by  means  of  an  appropri¬ 
ately  designed  linear  filter. 

Our  approach  to  the  interference  suppression  problem  has  been  greatly 
influenced  by  the  previous  work  of  Hsu  and  Giordano  [ 1 ] .  They  considered 
the  problem  of  narrow  band  interference  estimation  and  suppression  by 
means  of  two  linear  prediction  algorithms,  the  Burg  algorithm  |2,x),  and 
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the  Levinson  algorithm  |2,4|.  The  channel  through  which  tin'  PN  spread 
spectrum  signal  is  transmitted  was  assumed  to  he  nondispersi ve.  Results 
were  presented  on  the  effectiveness  of  the  linear  prediction  filter  in 
suppressing  the  interference.  Performance  was  measured  in  terms  of 
signal -to- noise  ratio  at  the  output  of  the  PN  correlator. 

Our  research  work  extends  the  results  obtained  by  llsu  and  Giordano 
on  filter  requirements  and  characteristics  in  single  and  multiple  fre¬ 
quency  band  interference.  In  addition  to  a  nondispersive  channel,  we 
consider  the  transmission  of  the  PN  spread  spectrum  signal  over  a 
channel  characterized  by  fading  and  multipath  (time  dispersion).  This 
serves  as  a  model  for  radio  channels  such  as  III-'.  The  existence  of  time 
dispersion  in  the  received  signal  necessitates  some  means  for  dealing 
with  this  type  of  distortion  at  the  receiver.  We  have  considered  the  use 

i 

of  an  adaptive  decision- feedback  equalizer  preceding  the  PN  correlator 
for  mitigating  the  effects  of  time  dispersion  due  to  multipath  and  the 
linear,  interference  suppression  filter. 

Section  II  of  this  report  presents  the  algorithms  for  estimating 
and  suppressing  narrow  band  interference  in  a  wideband  PN  spread  spectrum 
signal.  Performance  results  are  presented  in  Section  Til.  Before  con¬ 
cluding  this  section,  we  present  a  brief  description  of  the  mathematical 
model  of  the  PN  spioad  spectrum  binary  communication  system  which  is  used 

in  the  analysis  and  in  Monte  Carlo  simulation.  A  baseband  system 
is  used  throughout  this  report. 

1.1  Mathematical  Model  of  PN  Spread  Spectrum  Communication  System 

Transmitted  Signal 

The  transmitted  signal  is  generated  as  shown  in  the  block  diagram. 


B i na rv 
Data 


Modulator 


sit) 


PX  Code 
Cenerator 


The  number  of  P.\  chips  per  information  bit  is  I.,  ihus,  the  signal  for 
the  ktb  infonnation  bit  can  be  expressed  as 


bk(t)  =  i\j  l>(t  -  -'V 


(1.1) 


where  ip^  1  represent  the  output  sequence  from  the  P.\  code  generator  for 
the  kTb  information  bit  and  q(t)  is  a  rectangular  pulse  of  duration  t 
and  unit  energy.  Hie  total  transmitted  signal  may  be  expressed  in  the 
form 


Mtl  -  I  Ik  hk(t  -  kThl 


Cl  -  - 1 


where  l^i  represents  the  binary  information  sequence  and  =  l.rc  is  the 
bit  interval  (reciprocal  of  the  bit  rate). 

Cliannel 

Hie  fading  multipath  channel  is  modeled  as  having  discrete  multi- 
path  components  with  relative  delays  equal  to  multiples  of  the  chip 
duration  ;  .  The  impulse  response  may  be  expressed  as 

k 

hit;:)  -  j  <X  (  t  )  b  ( T  -  j  T  .  )  (1.0) 

;  - 1  1  *- 
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where  the  {cu(t)}  are  complex-valued,  statistically  independent,  narrow  band 
Gaussian  random  processes.  Since  the  multipath  components  cannot  be 
distinguished  with  a  resolution  better  than  x  ,  this  is  a  reasonable 
model  for  the  multipath.  The  fading  is  inherent  in  the  time  variations 
of  the  narrow  band  processes  (cij(t)}. 

The  signal  is  corrupted  by  additive  white  Gaussian  noise  n(t)  and 
by  narrow  band  interference  denoted  by  i(t).  The  narrow  band  inter¬ 
ference  is  modeled  as  consisting  of  either  a  number  of  CW  tones,  i.e., 

Q 

i(t)  =  l  Am  cos (2ir  ^t  +  <|>  )  ,  (1.4) 

m=l 

or  a  filtered  narrow  band  Gaussian  noise  process.  In  some  cases  it  may 
be  appropriate  to  model  the  interference  as  a  narrow  band  random  process 
that  arrives  at  the  receiver  through  another  statistically  independent 
fading  multipath  channel.  Although  we  have  not  considered  this  model  of 
the  interference  explicitly  in  our  analysis,  the  algorithms  presented 
in  Section  II  for  estimating  and  suppressing  the  interference  still  apply. 

Received  Signal 

From  the  description  given  above,  the  received  signal  has  the  form 
K 

r(t)  =  I  a.(t)  s(t  -  jx  )  +  i(t)  +  n(t)  (1.5) 

j  =  l  J 

The  receiver  attempts  to  suppress  the  interference  i(t)  and  then  to  re¬ 
cover  the  information  sequence  by  further  processing,  which  involves 
equalization  and  cross-correlation  with  a  replica  of  the  PN  sequence. 
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Estimation  and  Suppression  of  the  Interference 


The  estimation  and  suppression  of  the  narrow  band  interference  is 
accomplished  prior  to  signal  demodulation,  as  illustrated  in  general 
terms  in  the  block  diagram  shown  below. 


r(t) 


The  estimation  of  the  filter  coefficients  may  be  accomplished  by  means 
of  a  linear  prediction  algorithm  or  by  means  of  a  spectral  analysis  al¬ 
gorithm  based  on  the  Fast  Fourier  Transform  (FFT)  algorithm.  In  any 
case,  the  objective  is  to  design  an  adaptive  transversal  filter  that 
highly  attenuates  the  received  signal  in  those  frequency  bands  which  con¬ 
tain  strong  interference. 

Signal  Demodulator 

The  signal  at  the  output  of  the  interference  suppression  filter  is 
processed  by  an  adaptive  equalizer  operating  on  a  chip-by-chip  basis, 
followed  by  a  PN  correlator  which  employs  a  replica  of  the  transmitted 
PN  sequence.  Perfect  synchronization  of  the  PN  sequence  is  assumed. 

The  type  of  equalizer  employed  in  our  performance  evaluation  is  a 
decision- feedback  equalizer.  It  serves  to  reduce  the  detrimental  effects 
of  time  dispersion  due  to  the  interference  suppression  filter  and  the 
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In  this  section  we  present  a  number  of  algorithms  for  estimating  and 
suppressing  a  narrow  band  interference  embedded  in  a  wideband  PX  spread 
spectrum  signal.  The  algorithms  may  be  classified  into  two  general  cate¬ 
gories.  The  algorithms  in  the  first  category  employ  the  fast  Fourier 
Transform  (ITT)  algorithm  for  performing  a  spectral  analysis  from  which 
an  appropriate  transversal  filter  is  specified.  These  algorithms  may  be 
termed  nonparametric,  since  no  prior  knowledge  of  the  characteristics  of 
the  interference  is  assumed  in  forming  the  estimate.  The  algorithms  in 
the  second  category  are  based  on  linear  prediction  and  may  be  termed  para¬ 
metric.  That  is,  the  interference  is  modeled  as  having  been  generated  by 
passing  white  noise  through  an  all-pole  filter  f 2 1 . 

2.1  Interference  Suppression  Based  on  Nonparametric  Spectral  listimates 

'flic  basis  for  this  method  is  that  the  power  density  spectrum  of  the 
P.\  sequence  is  relatively  flat  while  the  spectrum  of  the  narrow  band 
interference  is  highly  peaked.  The  first  step  in  this  method  is  to 
estimate  the  power  spectral  density  of  the  received  signal.  The  spectral 
estimate  can  he  obtained  by  any  one  of  the  well-known  spectral  analysis 
techniques  described  in  1 5 J .  for  illustrative  purposes,  we  have  selected 
the  Welch  method  and  we  have  made  use  of  the  computer  program  listed  in 
[b|  to  generate  the  numerical  results  presented  in  Section  III. 

i ince  the  power  spectral  density  of  the  received  signal  is  estimated, 
the  interference  suppression  filter  can  be  designed.  A  transversal  filter 
is  an  appropriate  filter  structure  for  this  application,  since  we  desire 
to  use  a  filter  that  contains  ceros  in  the  frequency  range  occupied  by  the 


interference.  A  relatively  simple  method  for  designing  the  transversal 


filter  in  the  discretc-t ime  (sampled-data )  domain  is  to  select  its  dis¬ 
crete  Fourier  transform  (OFT)  to  be  the  reciprocal  of  the  square  root  of 
the  [lower  spectral  density  at  equally  spaced  frequencies.  To  elaborate, 
suppose  that  the  transversal  filter  has  K  taps.  The  problem  is  to  specify 
the  K  tap  coefficients  h(n)}  or,  equivalently,  the  HIT  llfk),  defined  as 

F- 1  - i  IT  nk 

llfk)  =  l  h (n)  e  ,  k  =  0,1 _ _ _  K- 1  (2.1) 

n=0 


The  I 'FT  1 1  ( k ) ,  k  =  0,1,  .  ..,  K-l  is  selected  as 


I  Ilk)  = 


•  CM  , 

■J  X  1 


e 


(2.2) 


where  i’(f),  0  <  f  <  R  ,  denotes  the  estimate  of  the  power  spectral 
density  and  R^.  denotes  the  sampling  rate,  which  is  normalized  to  unity. 

It  is  desirable  to  have  a  transversal  filter  that  has  linear  phase. 
This  can  be  achieved  if  the  impulse  response  h(n)  is  real  and  satisfies 
the  symmetry  condition 


h  ( n )  =  h  ( K  -  1  -  n ) 


(2.5) 


The  symmetry  condition  in  Fq.  (2.5)  is  satisfied  if  ll(k)  =  ll*(k  -  k)  . 
But  ll(k)  as  defined  in  Fq.  (2.2)  docs  satisfy  this  condition  since  P(f)  = 
P(RS  -  0 •  Hence,  h(n)  is  symmetric. 

In  effect,  the  filter  characteristic  obtained  from  Fq.  (2.2)  attempts 
to  approximate  an  inverse  filter  to  the  power  spectral  density.  That  is, 

-S 


the  interference  suppression  filter  attempts  to  whiten  the  spectrum  of 
the  incoming  signal.  Thiis,  the  filter  will  have  a  large  attenuation  in 
the  frequency  range  occupied  by  the  interference  and  a  relatively  small 
attenuation  elsewhere. 

There  is  an  alternative  filter  design  procedure  that  may  lead  to  a 
relatively  smaller  transversal  filter.  'Hie  method  simply  involves  the 
selection  of  the  position  of  the  zeros  so  as  to  obtain  an  appropriate 
set  of  notches  in  the  frequency  response  characteristic  of  the  filter. 

Ivc  have  not  investigated  any  ad  hoc  methods  for  selecting  the  zeros  of 
the  transfer  function,  since  the  linear  prediction  approach  described  in 
the  following  section  is  a  systematic  method  for  attaining  the  same  goal. 

2.2  Interference  Suppression  Ba_se_d_ _on_  Linear  Prediction 

In  contrast  to  the  nonparamet ric  spectral  analysis  method  described 
in  the  previous  section,  the  method  presented  in  this  section  for  esti¬ 
mating  the  narrow  band  interference  is  based  on  modeling  the  interference 
as  white  noise  passed  through  an  all-pole  filter.  That  is,  instead  of 
using  the  received  signal  to  estimate  the  spectrum  directly,  the  signal 
is  used  to  estimate  the  pole  positions.  This  estimation  is  accomplished 
by  means  of  linear  prediction.  An  estimate  of  the  power  spectral  density 
is  easily  obtained  from  the  all-pole  model.  However,  this  step  can  be 
omitted.  That  is,  the  power  density  spectrum  need  not  be  computed  ex¬ 
plicitly  for  the  purpose  of  designing  the  suppression  filter.  The  inter¬ 
ference  suppression  filter  is  simply  a  transversal  (all-zero)  filter 
having  zero  positions  that  coincide  with  the  estimated  pole  positions. 
Thus,  the  spectrum  of  the  signal  at  the  output  of  the  transversal  filter 


is  rendered  white. 


In  order  to  develop  the  mathematical  formulation  for  the  all-pole 
model,  we  assume  that  the  channel  through  which  the  signal  is  transmitted 
is  nondispersive  (no  multipath).  Let  sft)  denote  the  equivalent  lowpass 
transmitted  PN  spread  spectrum  signal  and  let  r(t)  denote  the  equivalent 
lowpass  received  signal.  The  latter  is  expressed  as 

r(t)  =  s(t)  +  i(t)  +  n(t)  (2.4) 

where  i(t)  denotes  the  narrow  hand  interference  and  n(t)  is  assumed  to  be 
a  sample  function  of  a  white  Gaussian  noise  process.  For  convenience,  we 
assume  that  r(t)  is  sampled  at  the  chip  rate  of  the  PN  sequence.  Thus, 
Lq.  (2.4)  can  be  expressed  as 


r(k)  =  s(k)  +  i(k)  +  n(k),  k  =  1,2,  ...  (2.5) 


U'c  assume  that  s(k),  i(k)  and  n(k)  are  mutually  uncorrelated. 

An  estimate  of  the  interference  i ( t )  is  formed  from  r(k).  Assume 
for  the  moment  that  the  statistics  of  i(t)  are  known  and  are  stationary. 
Then,  we  can  predict  i(k)  from  r(k  -  1),  r(k  -  2),  ...,  r(k  -  m) .  That 

i-% 


m 

i  ( k )  =  l  a,  r(k  -  i)  (2.6) 

?  =  1  4 


where  fa„I  arc  the  coefficients  of  the  linear  predictor.  It  should  be 
emphasised  that  Lq.  (2.6)  predicts  the  interference  hut  not  the  signal 
s  ( k ) ,  because  s(k)  is  uncorrclated  with  rik  -  ?.)  for  Z  =  1,  2,  ...,  m, 

'Tor  convenience,  our  treatment  oi  linear  prediction  is  based  on  real-valued 
signals.  The  extension  to  complex-valued  signals  is  straightforward. 
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as  a  consequence  of  the  sampling  being  done  at  the  chip  rate. 

The  coefficients  in  liq.  (2.0)  are  determined  by  minimizing  the  mean 
square  error  between  r(k)  and  ilk),  which  is  defined  as 

£(.m)  =  !'.|  r(k)  -  i(k)  j2 

=  l:|r(k)  -  l  a  r(  k  -  l)  |“  (2.7) 

Jl*l 

Minimization  of  E  with  respect  to  the  predictor  coefficients  {a  ')  can 
he  easily  accomplished  by  invoking  the  orthogonality  principle  in  mean 
square  estimation  f 7 1 .  This  leads  to  the  set  of  linear  equations 

m 

I  a  p(k  -  £)  =  p(k)  ,  k  =  1,2,  ...,  m  (2.8) 

n=i 

where 

P(k )  =  bfr(m)  r(k  +  m) )  (2.9) 

is  the  autocorrelation  function  of  the  received  signal  r(k)  .  The 
equations  in  (2.8)  are  usually  called  the  Yule-Walker  equations  [2]. 

They  can  be  solved  efficiently  by  means  of  the  Levinson  algorithm  f 4 , 2 1  . 

The  Levinson  algorithm  is  an  order- recursive  method  for  solving  liq. 
(2.8).  That  is,  it  solves  for  the  coefficients  of  an  m- order  predictor 
recursively  from  the  coefficients  of  an  (m  -  1) -order  predictor.  Starting 
with  a  first-order  predictor  and  introducing  another  subscript  in  the  pre¬ 
diction  coefficients  to  indicate  the  order,  we  have 
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and,  hence, 


Pin)  an  =  Pin 


*11 


pm 

pfoT 


(2.10) 


For  a  second-order  predictor,  the  two  equations  obtained  from  Hq.  (2.8) 
are 


p(0)  a,1  +  p(l)  a??  =  p(l) 

Pfl)  a21  +  P CO)  a22  =  0(2)  (2.11) 

The  first  equation  in  (2.11)  can  he  used  to  solve  for  a 2j.  By  substi¬ 
tuting  a^  for  p(l)/DfO)  we  obtain 

^■>1  =  a^2  -  a71  (2.12) 

Thus,  a,j  is  related  to  a^.  \'ext,  a^  can  be  solved  from  the  second 
equation  in  (2.11)-  Bv  using  Tq.  (2.12)  to  eliminate  a^  from  the  second 
equation,  we  obtain  a77  in  the  form 

P (2 )  -  p(l)  an 

a22  =  p(0)  -  r(irajY  (2,1^ 

Therefore,  Hq.  (2.13)  is  used  to  solve  for  a,,  and,  then,  a-,.j  is  ob¬ 
tained  from  Hq .  (2.12).  This  is  the  Levinson  recursion  for  order  two. 

In  general,  it  can  be  shown  [2,4 1,  that  t  lie  Hcvinson  recursion  for  the 
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coefficients  of  the  m- order  predictor  are 


amk  am- 1  k  amm  am-l  m-k 


,  k  =  1,2,  . . . ,  m- 1 


14) 


where 


pfm)  -  f 


-m-1  -in-1 


min  Pin  -  f  ,  ar  ; 

-ill- 1  -m- 1 


(2.15) 


The  vectors  h  ,  a  and  ar  in  Tq.  (2.15)  are  defined  as 
-n  -n  1 


"  *  in)  " 

~  a  ”1 

nn 

f  a  .  “I 
nl 

= 

P(n-l) 

r  _ 

’  ^-n 

a  i 

nn  - 1 

’  an 

a  _ 
n2 

_  Pill 

-  an  1  - 

a 

L  nn  -1 

(2.16) 


r 

IVc  note  from  f.qs .  (2.16)  that  the  vector  a^  is  simply  the  vector  a 
in  reverse  order,  furthermore,  if  we  express  Hq.  (2.8)  in  the  matrix 
form 


where  1^  is  the  (tnxm)  autocorrelation  matrix,  the  vector  in  liqs.  (2.16) 
is  simply  the  vector  h  in  reverse  order.  These  relations  will  he  used 
later  in  our  discussion  of  the  lattice  realization  of  the  prediction 


fi Iter. 

The  minimum  mean  square  error  is  a  measure  of  the  effectiveness  of 


the  prediction  filter.  The  expression  for  £'^^(111)  is  easily  shown  to 
be 


m 


J'«)  =  ”(0)  -  I  amk  f  <k) 


m  l  n 


k=l 


(2.18) 


-\  recursion  relation  can  also  be  obtained  for  £ u^n (m) .  Using  hq.  [2.18) 
we  have 


m- 1 

;Jm)  =  pfO)  =  I  amkp(k)  -  a  o(m) 
k=l 


nun 


(2.19) 


If  we  substitute  for  a  ,  in  F.q.  (2.19)  the  Levinson  recursion  given  in 

mk 

Hq.  (2.14),  and  rearrange  the  terms,  we  obtain 


m- 1  m 

W>  ■  toco)  -  ^  vik  Vi  m-k  »»>: 


nun 


(m  '  n  -  amm[p(-m)  '  ^-1  (2'20) 


-m- 1  --m- 1 


From  Fq.  (2.15),  we  note  that 


-  Pim)  ~  ^-1  Vi 

'mm  "  A'minCm  -  D 


(2.21) 


Hence,  Hq.  (2.20)  becomes 


F  .  (m)  =  (1  -  a"  )  E  .  (m  -  1) 
m  1  n  iron  m  1  n 


!.  22) 
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The  result  in  hq.  (2.22)  implies  that  la  1. 

Once  the  prediction  coefficients  arc  determined,  the  estimate  i ( k ) 
of  the  interference,  given  hv  hq.  (2.(i),  is  subtracted  from  r ( k J  and 
the  difference  signal  is  processed  further  in  order  to  extract  the  digital 
information.  Thus,  the  equivalent  transversal  filter  for  suppressing  the 
interference  is  described  by  the  transfer  function 


Vs 


=  1 


L 

k=l 


mk 


(  2 . 2  a ) 


where  :  '  denotes  a  unit  of  delay.  The  corresponding  all -pole  model  for 

the  interference  signal  is  1/A  (z). 

m 

The  solution  of  hq.  (2.K)  for  the  coefficients  {a  ,  }  of  the  pre¬ 
diction  filter  requires  knowledge  of  the  autocorrelation  function  r(k). 

In  practice,  the  autocorrelation  function  of  ilk)  and,  hence,  r(k)  is 
unknown  and  it  may  also  be  slowly  varying  in  time.  Consequently,  one 
must  consider  methods  for  obtaining  the  predictor  coefficients  directly 
from  the  received  signal  i r ( k ) } .  This  may  be  accomplished  in  a  number 
of  ways.  In  this  investigation,  three  different  methods  were  considered. 
In  all  cases,  we  obtained  the  predictor  coefficients  by  using  a  block  of 
N  samples  of  frlklk  Hie  three  methods  are  described  in  the  following 
sect  ion . 


2. A  Algorithms  for  Computing  the  Prediction  Coefficients 

Since  the  autocorrelation  function  o(k)  is  not  known  r  vtrvi ,  we 
shall  describe  three  algorithms  for  computing  the  prediction  coefficients 
from  the  received  signal  rfk).  We  assume  throughout  that  a  block  of  N 
samples  of  r(k)  are  available. 
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Direct  Application  1  of  the  Levinson  Algorithm 

The  first  method  is  simply  based  on  the  direct  estimation  of  pfkj 
from  t he  block  of  X  samples.  The  estimate  of  rfk)  is 

X-k 

P ( k)  =  l  r (n)  r(n  +  k)  ,  k  =  0,1,  ...,  m  (2.24) 

n=  1 

fhe  estimate  p ( k)  may  then  be  substituted  in  Lq.  (2.8)  in  place  of  o(k) 
and  the  Levinson  algorithm  can  be  used  to  solve  the  equations  efficiently. 
Thus,  the  recursive  relations  given  in  the  previous  section  apply  with 
p(k)  replaced  by  p ( k ) . 

Burg  Algorithm 

The  second  method  considered  for  obtaining  the  prediction  coefficients 
is  the  Burg  algorithm  [2,31  -  Basically,  the  Burg  algorithm  may  be  viewed 
as  an  order- recursive  least  squares  algorithm  in  which  the  Levinson  re¬ 
cursion  is  used  in  each  iteration.  To  be  specific,  we  begin  with  the 
determination  of  the  coefficient  in  a  first-order  predictor,  based  on  the 
method  of  least  squares. 

The  performance  index  used  by  Burg  is  the  mean  (time-average)  square 
error,  which,  for  a  first-order  predictor,  is  defined  as 

\' 

f  Bf  1 )  =  ),  1 1  r  (  i )  -  a-jj  r(i  -  1)]“  +  [r(i  -  1)  -  a]  1r  f  i)  ]  “ } 

i=2 

X 

=  l  |  f j  (  i )  +  b“(  i )  |  (.2.25) 

i  =  2 

file  first  term  in  t fie  sum  represents  t lie  error  in  a  first-order  predictor 
operating  on  the  received  signal  in  the  forward  direction  and,  hence,  it 


is  called  the  forward  error  f^(i).  The  second  term  represents  the  error 
in  a  first-order  predictor  operating  on  the  data  in  the  reverse  direction 
and,  hence,  it  is  called  the  backward  error  ( i ) .  Minimisation  of 
fp[l)  with  respect  to  a^  yields  the  predictor  coefficient  n^  in  the 
fonn 

X 

2  l  r(i)  r ( i  -  1) 

a  = _ Azi: _ _ _  (->  '>(,1 

all  N  7  ,  U  1 

l  I r“ ( i )  +  r"Ci  -  D I 

i=2 

for  the  second-order  predictor  we  have  the  performance  index 
X  t 

£B(2)  =  l  lf;(i)  +  b;( i ) I  (2.27) 

i=5 

where  f?(i)  and  h7(i)  denote  the  forward  and  backward  errors,  respectively, 
which  are  defined  as 

2 

f i (i )  =  r(i)  -  l  a7,  r(i  -  k) 

"  k=l  ZK 

7 

b7(i)  =  r(i  -  2)  -  l  a,.  r(i  -  2  +  k)  (2.28) 

k=l  “ 

We  use  the  Levinson  recursion  given  in  Lq.  (2.12)  to  express  f,(i)  and 
h,(ij  in  tenns  of  a^  and  a-,-,,  where  a^  is  the  coefficient  of  the  first- 
order  predictor,  lor  f/i)  we  obtain  the  relation 


f2(.i)  =  r(i)  -  a}1  «'(i  *  D  -  a?2|r(i  -  2)  -  r(i  -  1J] 


Ihc  general  Levinson  recursion,  given  in  Lq .  (2.11),  is  used  in  liqs.  (2.5a) 
to  express  f  (i)  and  b  (;)  in  terms  of  the  predictor  coefficients  at  the 
(in  -  1)  iteration.  The  result  of  this  substitution  is 


1' 

in 


(i) 


f  .  (  i ) 
m- 1 


nun  in- 1 


(i 


1) 


h  (i) 


m- 1 


fi  -  1) 


mm 


fm-l(i) 


(2.34) 


Upon  substituting  liqs.  (2.34)  into  liq.  (2.32J  and  minimizing  f„(m)  with 

respect  to  a  we  obtain 
mm 


a 

mm 


N 

2  )  f  .  ( i )  b  .  ( i  -  1 ) 

»-i  _ 

.  I  |fi. ini  *  -  n 

i=m+ 1 


(2.35) 


A  recursion  relation  can  also  be  obtained  for  the  minimum  mean 
square  error.  By  substituting  liqs.  (2.34)  into  liq.  (2.32),  expanding  the 
squared  terms  and  using  liq.  (2.35),  we  obtain  the  recursive  relation 


’  .  (in)  =  (1  -  az  )  F,,  .  (m 

B  nun  mm  B  nun 


1) 


(2.3  b) 


The  recursive  relations  in  I  qs .  (2.34),  (2.35)  and  (2.3b)  constitute 
the  Burg  algorithm  for  computing  the  predictor  coefficients  and  the 
minimum  mean  square  error.  The  estimate  of  the  power  spectral  density 
obtained  by  means  of  the  Burg  algoritlun  using  a  predictor  of  order  m  is 
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pm 


1 


L\,  .  (m) 

B  in  in  ’ 

m 

I  a  i  o  -1 
k=l  mk 


(2.37) 


[.cast  Squares  Algorithm 

As  we  indicated  in  the  discussion  above,  the  Burg  algorithm  is 
basically  a  least  squares  algorithm  with  the  added  constraint  that  the 
predictor  coefficients  satisfy  the  Levinson  recursion.  As  a  result  of 
this  constraint,  an  increase  in  the  order  of  the  predictor  requires 
only  a  single  parameter  optimization  at  each  stage.  In  contrast  to  this 
approach,  we  shall  now  describe  an  unconstrained  least  squares  algorithm. 
That  is,  the  algorithm  computes  the  optimum  predictor  coefficients,  in 
the  sense  of  least  squares,  at  each  stage  of  the  iteration. 

As  in  the  Burg  algorithm,  we  minimize  the  mean  square  error  perform¬ 
ance  given  in  bq.  (2.32).  The  forward  and  backward  prediction  errors 
arc  defined  in  hqs.  (2.33).  The  global  minimization  of  bq.  (2.32)  with 
respect  to  the  set  of  predictor  coefficients  (a  jJ  yields  the  set  of 
linear  equations 

m 

l  a  ,  f>U,k)  =  <PU,0)  ,  e  =  1,2,  ...,  m  (2.38) 

k=  1  mK 


where 


N 

<+>(£., k )  =  l  |  r(  i  -  k )  r(  i  -  i)  +  r ( i  -  m  +  k)  r(i  -  m  +  £)  ] 
i  =m+ 1 

(2.39) 


The  linear  equations  in  (2.38)  can  lie  expressed  in  matrix  form  as 
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where  the  matrix  0  is  an  (mxni)  autocorrelation  matrix  with  elements 

nil 

J  ■'*  { k ,  i )  !  and  is  an  m-dimensional  vector  with  elements  4(?,,(l), 
k  =  1,2,  ....  in.  'Hie  matrix  v  is  symmetric.  However,  in  contrast  to 
the  autocorrelation  matrix  II  in  I iq .  (2.17),  which  is  Toepliti,  the 
matrix  is  not  Toeplitc.  Consequently,  the  Levinson  algorithm  cannot 
be  used  to  solve  Hq.  (2.40)  recursively.  In  spite  of  the  fact  that  4 
is  not  Toeplit:,  it  is  still  possible  to  derive  a  recursive  algorithm 
for  the  predictor  coefficients  based  on  the  least  squares  performance 
index.  Morf  et  al  [8-11]  have  developed  such  a  recursive  least  squares 
algorithm  that  not  only  allows  one  to  recursively  increase  the  order  of 
the  predictor,  but  also  allows  one  to  update  the  predictor  coefficients 
recursively  in  time  for  a  predictor  of  any  given  order.  A  detailed 
development  of  this  recursive  least  squares  algorithm  lias  been  given  by 
Pack  and  Sa tori us  [12]  and,  for  the  sake  of  brevity,  will  not  be  repeated 
here.  An  order- recursive  version  of  this  algorithm  has  also  been  de¬ 
scribed  recently  by  Marple  [13].  The  point  that  we  wish  to  make  is  that 
the  prediction  coefficients  based  on  the  least  squares  criterion  can  be 
solved  efficiently  by  means  of  an  algorithm  that  is  both  recursive  in 
order  and  in  time.  The  order- recursive  part  of  the  algorithm  fits  the 
lattice  formulation  described  briefly  in  the  following  section. 

To  conclude  this  section,  we  observe  that  the  minimum  mean  square 
error  in  t he  least  squares  solution  can  be  expressed  as 


(2.41) 


N' 

E.  c  .  (m)  =  5'  [f  (i)  r(i)  +  h  (i)  rfi  -  m)l 

I.S  nun  .  L  ,  L  m  in  1 

i=m+l 

in 

=  <{.(0,0)  -  l  a  .  <Mk,d| 

k=l 

As  a  final  comment  we  mention  the  well-known  property  of  the  unconstrained 

least  squares  solution,  namely,  that  the  resulting  prediction  coefficients 

do  not  necessarily  vield  a  minimum  phase  filter  A  ( z) .  In  our  case  and  in 
■  ■  m 

some  other  practical  applications  of  linear  prediction  this  is  not  a 
problem. 

2 . 4  Lattice  Structure  for  Linear  Prediction  filters 

In  this  section  we  demonstrate  that  the  transversal  filter  with 
transfer  function 


A  (z)  =  1 
nr 


m 


l 

k=l 


(2.42) 


where  the  (a  are  the  prediction  coefficients,  can  also  be  realized  as 
a  lattice  filter.  As  we  shall  observe  below,  the  lattice  filter  structure 
has  a  number  of  properties  that  may  prove  desirable  in  a  practical  imple¬ 
mentation  of  the  interference  suppression  filter. 

The  starting  point  for  this  development  is  the  Levinson  recursive 
algorithm  for  the  predictor  coefficients  {a^.}  given  in  liq.  (2.14).  If 
we  substitute  this  relation  into  liq .  (2.42)  and  rearrange  terms,  we 
obtain 
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Thus,  I’q.  (2.4b)  gives  the  recursive  relation  for  Gm ( 3 ) . 

Let  fin  and  b  (i)  denote  the  output  sequences  of  the  filters  A  (2) 

and  (I^fc),  respectively,  for  an  input  sequence  x(i).  Since,  I-'m( z)  = 

A  (2)  X(z)  and  B  (2)  =  Cl  (z)  X  (2)  ,  it  follows  from  T.qs.  (2.43)  and 

(2.45)  that  1  (2)  and  B  (2)  satisfy  the  recursive  relations 
’  m  J  m  J 
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where  initially  we  have  M  (2)  =  Bq(z)  =  X(z).  Alternatively,  in  the  time 
domain  the  relations  in  Hqs.  (2.47)  arc 
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with  initial  conditions  f  (i)  =  bp(i)  =  x(i). 

flic  recursive  relations  in  f.qs.  (2.48)  describe  a  lattice  filter  as 
shown  in  figure  2.1.  In  our  case  the  input  to  the  lattice  is  the  re¬ 
ceived  signal  sequence  r(i)  and  the  desired  output  is  f  (i).  Mach  stage 
of  the  lattice  has  two  inputs  and  two  outputs.  The  two  inputs  and  out¬ 
puts  are  the  forward  ;md  backward  prediction  errors  defined  previously. 
These  errors,  usually  called  the  residuals,  satisfy  a  number  of  inter¬ 


esting  properties  [11],  which  we  state  here  without  proof. 

first  of  all,  there  is  the  orthogonality  property  between  the 


residuals  ;md  the  input: 


h  [  i'  (  i  J  r  (  i  -  k  J  J  =  0  ,  1  <  k  <  m 

1  m  1  ’  - 


li|b  (i)  r(  i  -  k)  ]  =  0  ,  1  <_  k  £  m-1 


(2.49) 


These  relations  are  simply  a  restatement  of  the  orthogonality  principle 
in  linear  me  cm  square  estimation  ( 7  J .  Secondly,  the  backward  residuals 
arc  self-orthogonal,  in  the  sense  that 


H[h  ( i)  b  (i) ]  =  e  ■  (m)  c 

1  m  J  n  J  nun  '  mn 


(2. SO) 


where  F.  ■  (in)  is  the  minimum  mean  square  error  for  an  m-stage  lattice  or 
min  ’ 

predictor,  and  5  is  the  Kronecker  delta.  This  orthogonal itv  relation 
means  that  successive  stages  of  the  lattice  are  decoupled  statistically. 
In  terms  of  the  forward  residuals,  the  expression  for  the  minimum  mean 
square  error  is 


E  ■  (ill)  =  i:[f  (i) 
min  1  m J 


(2.51) 


Two  additional  expressions  for  !■'  ■  (m)  arc 
1  nun 


■'  .  (m)  =  l.|  f  (  i )  r(  i ) 
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Finally,  the  prediction  coefficients  {a^}  computed  at  each  stage  of  the 
lattice  are  related  to  the  cross-correlation  between  the  forward  and 
backward  residuals.  Specifically,  we  have 
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In  using  the  lattice  structure  instead  of  the  transversal  structure 
for  the  interference  suppression  filter,  one  must  compute  the  lattice 
gains  {a,jJ  as  specified  by  liq.  (2.15)  or  liq.  (2.21),  in  conjunction  with 
the  Levinson  recursion  (2.14).  In  this  computation,  the  estimate  (time- 
average)  of  the  autocorrelation  function  p(k)  is  used  in  place  of  the 
statistical  average  p(k). 

The  Burg  algorithm  as  described  by  the  relations  in  lips.  (2.54), 
(2.55)  and  (2.5b)  is  basically  a  lattice  implementation  of  the  linear 
predictor.  The  lattice  stages  are  specified  by  the  recursive  relations 
in  litis .  (2.54)  for  f  ( i  )  and  b  (i),  with  the  lattice  gains  {a  }  given 
by  l.q.  (2.55).  It  is  interesting  to  note  that  liq.  (2.55),  which  was  ob¬ 

tained  by  performing  a  least  squares  optimization,  is  the  time-average 
equivalent  of  the  statistical  average  given  by  lie).  (2.54).  In  fact,  the 
forward  and  backward  residuals  have  identical  statistical  mean  square 
values.  Consequently,  the  denominator  in  liq.  (2.55)  is  equivalent  to 
twice  the  statistical  moan  square  value.  Thus,  the  factors  of  two  in  the 
numerator  ;uul  denominator  cancel . 

from  a  computational  viewpoint,  the  Burg  algorithm  has  the  advantage 


that  the  lattice  gains  arc  computed  directly  from  the  forward  and  back¬ 
ward  residuals,  as  indicated  by  the  relation  in  liq.  (2.55),  whereas  in 


the  direct  application  of  the  Levinson  algorithm  one  must  first  compute 
the  estimate  of  the  autocorrelation  function  from  the  data.  As  a  con¬ 
sequence,  the  Burg  algorithm  is  computationally  more  efficient. 

The  least  squares  algorithm  described  in  the  previous  section  can 
also  be  formulated  in  terms  of  a  lattice  filter  structure.  Such  a 
formulation  was  developed  by  Morf  et  al  [8-111.  for  a  clear  tutorial 
presentation  of  the  lattice,  recursive  least  squares  formulation  the 
interested  reader  is  referred  to  the  report  by  Pack  and  Satorius  [12]. 

2.5  (leneralizat ion  of  Linear  Prediction  Algorithms 

Up  to  this  point,  our  treatment  of  linear  prediction  algorithms 
has  been  limited  to  one-step  prediction.  In  the  PN  spread  spectrum 
problem  under  consideration,  one-step  linear  prediction  is  appropriate 
if  the  received  signal  is  sampled  at  a  rate  of  one  sample  per  PN  chip. 

In  such  a  case  the  predictor  forms  an  estimate  of  the  narrow  band  inter¬ 
ference  and  is  insensitive  to  the  presence  of  the  desired  signal.  As  a 
consequence,  the  desired  signal  is  not  suppressed. 

If,  for  some  reason,  the  received  signal  is  sampled  at  a  higher  rate, 
say  M  samples  per  chip,  then  successive  s;uiiples  of  the  desired  signal  arc 
highly  correlated.  In  this  case,  the  one-step  linear  predictor  will 
attempt  to  predict  and  suppress  not  only  the  narrow  band  interference 
but  also  the  desired  signal.  Since  this  situation  is  undesirable,  the 
remedy  is  to  employ  an  M-step  linear  predictor,  which  can  be  implemented 
as  illustrated  in  figure  2.2. 

A  similar  problem  arises  in  attempting  to  suppress  a  narrow  band 
interference  in  t he  presence  of  resolvable  signal  multipath.  If  the 


multipath  components  in  the  l’.\  spread  speetjum  signal  span  a  time  inter¬ 
val  of  M  samples,  a  one-step  linear  predictor  of  order  M  or  greater  will 
attempt  to  estimate  ;uul  suppress  the  multipath  components  in  the  received 
signal.  Since  the  multipath  signal  components  represent  a  form  of  signal 
diversity,  their  suppression  by  the  linear  predictor  may  be  undesirable. 

To  avoid  this  situation  an  M-step  linear  predictor  can  be  employed.  In 
the  case  of  multipath,  however,  the  major  problem  with  this  approach  is 
that  if  M  is  large  the  interference  components  i(k)  and  i (k  -  M  -  n) , 
n  =  0,1,  ...,  m-1  arc  highly  dccorrelatcd.  As  a  consequence,  th  •  estimate 
of  the  interference  is  poor  and  so  is  the  performance  of  the  suppression 
filter.  When  the  interference  is  much  stronger  than  the  desired  signal, 
the  multipath  components  of  the  signal  are  completely  masked  by  the 
interference.  In  this  case  a  one-step  predictor  is  appropriate,  since 
the  predictor  responds  to  the  interference  and  is  relatively  unaffected 
by  the  much  weaker  multipath  components  of  the  signal.  This  situation  is 
demonstrated  in  Section  III  by  some  numerical  results. 

In  spite  of  this  one  shortcoming  of  the  M-step  predictor,  it  may 
still  be  appropriate  to  use  it  in  some  situations,  as  for  example ,  in 
the  case  of  multiple  samples  per  chin  for  this  reason  we  briefly  outline 
the  appropriate  algor  it  Inns  for  performing  M-step  prediction  and  inter¬ 
ference  suppression . 

first,  we  adopt  a  statistical  approach.  The  M-step  linear  predictor 
of  order  m  is 

i’i 

id)--  l  c,  rit  -  k  -  M  +  1)  (2 . 55) 

k=l 

Tu 


and  the  corresponding  mean  square  error  is 


A' (in )  =  Id  |  r  (t )  -  i  ( t  J  |~]  1 2 . 5< >  J 

Minimization  of'  /■ 1  ( in )  with  respect  to  the  predictor  coefficients  {c.  } 
yields  the  set  of  linear  equations 

m 

l  Ci.  p  ( x  -  k)  =.(?,  +  M  -  1)  ,  1  =  1,:,...,  m  (2.57) 

k=l  K 

The  difference  between  this  set  of  linear  equations  and  the  set  given 
in  Lq.  (2.8)  for  the  one-step  predictor  are  the  terms  in  the  right-hand 
side  of  the  equations.  The  solution  of  Hq.  (2.57)  requires  a  general¬ 
ization  of  the  Levinson  algorithm  which  incorporates  the  recursion  in 
Lq.  (2.14)  and  includes  a  second  recursive  relation.  It  is  straight¬ 
forward  to  show  that  the  second  recursive  relation  is 
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k  =  1,2,  ....  m-1  (2.58) 


where  the  coefficients  c  is  given  bv  the  expression 
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The  vectors  a^  and  in  l.q.  (2.5‘d)  were  previously  defined  in  l.qs.  (2.1b), 

and  the  vector  c  is  the  n -dimens ional  vector  of  the  coefficients  {c  ,  1 , 
n  nk 

k  =  1,2,  ...,  n.  The  coefficients  (a  ;■  are  given  by  Lq.  (2.15).  1  inn  1 1  y , 
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the  initial  conditions  arc 


Substituting  Hqs.  (2.(i(>)  into  liq.  (2.t>5)  and  minimizing  the  resulting 


\ ... ( m )  with  respect  to  c  vields  t he  equation  for  c  in  the  form 
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Thus,  hqs.  (2. (>(>),  IJ.67)  along  with  (2.54)  and  (2.35)  constitute  a 
generalized  version  of  the  Burg  algorithm  appropriate  for  M-step  prediction, 
figure  2.4  illustrates  a  single  stage  in  a  lattice  implementation  of  this 
algorithm. 

If"  we  drop  the  constraint  that  the  predictor  coefficients  fc^ }  satis¬ 
fy  the  generalized  Levinson  algorithm,  and  simply  minimize  the  mean  square 
error  in  f.q.  (2.(i5)  with  respect  to  the  entire  set  of  coefficients,  we 
obtain  the  solution  for  the  least  squares  M-step  predictor.  The  appropri¬ 
ate  set  of  equations  for  the  coefficients  arc 


m 

}'  C  i,  41  f , k )  =  <p(*,0l  ,  f  =  1,2,  ...,  m  (2.bS) 
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where,  by  definition. 
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The  minimum  mean  square  error  is 


The  solution  to  the  above  equation  can  he  formulated  in  terms  of 
recursive  relations,  which  are  simply  generalizations  of  the  work  of 
Mori'  et  al  1 8 - 1 1 1  and  Pack  and  Satorius  [  1 2 J  .  The  recursive  relations 
represent  a  lattice  formulation  of  the  least  squares,  M-step  prediction 
problem . 

When  the  statistical  characteristics  of  the  received  signal  sequence 
r(i)  are  changing  with  time  (nonstationary  signal)  great  care  must  be 
taken  in  selecting  the  block  size  of  the  data  record  from  which  the  pre¬ 
diction  coefficients  are  determined.  In  particular,  N  must  be  small 
relative  to  the  number  of  samples  over  which  a  significant  change  occurs 
in  the  statistics  of  the  signal.  This  is  a  situation  where  the  parametric 
methods  based  on  linear  prediction  show  their  superiority  over  the  con¬ 
ventional  nonparamet ric  spectral  estimation  methods.  That  is,  for  the 
s.’irne  spectral  resolution,  the  parametric  methods  require  a  much  smaller 
set  of  samples  in  comparison  with  conventional  spectral  estimation.  Con¬ 
sequently,  the  parametric  method  will  accommodate  more  rapidly  changing 
signal  statistics. 

A  convenient  method  for  dealing  with  the  slowly  time-variant  proper¬ 
ties  of  the  received  signal  is  to  emplov  a  weighting  function  in  the 
performance  index,  which  places  more  emphasis  on  errors  from  current 
data  and  less  on  errors  due  to  past  data.  A  simple  weighting  function 
for  accomplishing  this  type  o!  weighting  is  an  exponential. 


In  the  context  of  one-step  prediction,  Morf  et  al  [8-11 1  and  Pack 


and  Satorius  [ 1 2 |  have  presented  a  recursive  least  squares  algorithm 
for  an  exponentially  weighted  mean  square  error.  Also  in  the  context  of 
one-step  prediction.  N'uttall  |18|  has  generalised  the  Burg  algorithm  by 
including  a  weighted  mean  square  error  performance  index. 

This  concludes  our  discussion  of  algorithms  for  narrow  hand  inter¬ 
ference  estimation  and  suppression.  In  the  next  section  we  present  a 
number  of  numerical  results  that  illustrate  the  effectiveness  of  inter¬ 
ference  suppression. 
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This  section  deals  with  the  performance  of  the  interference  sup¬ 
pression  filter.  First,  we  determine  the  improvement  provided  by  inter¬ 
ference  suppression  as  measured  in  terms  of  the  signal -to-noise  ratio 
(SNR)  at  the  output  of  the  !’.\  correlator.  The  output  SNR  is  by  far  the 
most  convenient  performance  index  for  obtaining  numerical  results.  This 
performance  index  is  used  to  assess  the  improvement  in  performance  ob¬ 
tained  by  an  interference  suppression  filter.  In  addition,  we  describe 
a  number  of  other  characteristics  of  the  interference  suppression  filter, 
including  its  frequency  response  and  the  location  of  its  zeros.  Finally, 
we  present  some  Monte  Carlo  simulation  results  on  the  perfonnance  of 
the  receiver  as  measured  in  terms  of  the  probability  of  error.  A  two- 
path,  Rayleigh  fading  channel  is  used  in  the  simulation.  The  receiver 
consists  of  an  interference  suppression  filter  followed  by  a  decision- 
feedback  equalizer  and  a  I’N  correlator. 

3.1  SNR  Improvement  Factor  Resulting,  from  Interference  Suppression 

In  order  to  demonstrate  the  effectiveness  of  the  interference  sup¬ 
pression  algorithms,  we  shall  compare  the  performance  of  the  receiver 
with  and  without  the  suppression  filter.  Since  the  channel  characteristic 
is  not  an  issue  in  this  type  of  comparison,  we  assume  that  the  channel  is 
ideal,  i.e.,  nondi spersive .  Consequently,  the  received  signal,  sampled 
at  the  chip  rate,  can  be  represented  as 

r(k  t  =  p(k  )  *■  i  ( k  !  +  n(  k )  ,  k  =  1,2,  ...  (3.1) 
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where  the  binary  sequence  { p ( k ) }  represents  the  P\  chips,  f i f k  J  >  repre¬ 
sents  the  sequence  of  samples  of  the  narrow  band  interference  and 


{n(k|}  represents  the  sequence  of  wideband  noise  samples. 

Let  the  impulse  response  of  the  interference  suppression  filter  be 
denoted  by  { h ( k ) } ,  k  =  0,1,  K.  In  the  case  of  the  one-step  pre¬ 

dictor  of  order  m,  wc  have  h(0)  =  1,  h(k)  =  -a^,  and  K  =  m.  For  an 
M-step  predictor  (M  >  1),  we  have  h(0)  =  1,  h(k)  =  0  for  1  <  k  <  M-l, 
hfk  +  M  -  1)  =  -c  .  for  k  =  1,2,  ...,  m  and  K  =  m  +  M  -  1.  The  input  to 
the  filter  is  r(k)  and  its  output  is 

k 

yfk)  =  l  hm  r(k  -  £)  ,  k  =  1,2,  ... 

k'.=0 

k 

=  l  hm|p(k  -  i)  +  i(k  -  Z)  +  n(k  -  &)|  f 3 . 2 ) 

f..=o 

It  is  apparent  from  Fq.  (3.2)  that  the  time-dispersive  characteristic  of 
the  interference  suppression  filter  results  in  inter-chip  interference 
which  can  be  mitigated  by  use  of  some  form  of  equalization  (16] .  In  the 
following  computation,  however,  no  equalizer  is  employed.  Instead,  the 
output  of  the  interference  suppression  filter  is  fed  directly  to  the  PN 
correlator.  Thus,  the  output  of  the  PN  correlator,  which  is  the  decision 
variable  for  recovering  the  binary  information,  is  expressed  as 

!. 

U  3  i  v( k )  p(k)  (3.31 

k  =  1 

where  I.  represents  the  number  of  chips  per  information  bit  or  the 


-40- 


processing  gain. 


substituting  liq. 

(3.2)  into  liq.  (5.3] 

i  ,  the 

decision  variable 

xpressed  in  the 

form 
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i  k 
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K  L  L  K 
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L  K 

+  [  I  hrn  p(k)  n(k  -  q  (3.4) 

k  =  l  Ji=o 

The  first  term  in  the  right-hand  side  of  liq.  (3.4)  represents  the  desired 
signal  component;  the  second  term  represents  the  self-noise  caused  hv 
tie  dispersive  characteristic  of  the  filter;  the  third  term  represents 
the  residual  narrow  band  interference  at  the  output  of  the  PN  correlator 
and  the  last  term  represents  the  additive  wideband  noise. 

For  the  comparison  that  we  wish  to  make,  the  S.\R  at  the  output  of 
the  l’\  correlator  is  a  mathematically  tractable  performance  index.  To 
determine  the  expression  for  the  SNR  we  must  compute  the  mean  and  variance 


of  IJ.  We  assume  that  the  binary  1\\  sequence  is  white,  the  interference 

i ( k )  has  zero  mean  and  autocorrelation  function  p.(k),  and  the  additive 

noise  n(k)  is  white  with  variance  Then,  the  mean  of  U  is 

n 

f.(  Ul  =  1.  (3.5) 


and  the  variance  is 
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K  7  K  k 

var(U)  =  1.  >  h"(?  )  +  I.  J  }  h f  •  )  h(m)  p.  (£  -  m) 
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I  he  first  tern  on  the  ri.uht-hand  side  of  the  expression  for  the  variance 
represents  the  mean  square  value  of  the  self-noise  due  to  the  time  dis- 
pet-ien  introduced  by  the  interference  suppression  filter.  The  second 
ten:  is  t hi  mean  square  value  el  the  residual  narrow  band  interference. 

he  last  t(  rn  i tlie  me-m  squat  e  value  of  the  wideband  noise. 

file  S\h  it  the  output  of  the  correlator  is  defined  as  the  rat  i  o  of 
the  square  ot  the  mean  to  t  lie  variance.  Thus, 

1. 

s  ,,  p  K  k  k 

hi-1  *  hi  -  i  him)  .  .  (  -  m)  +  J  h"(c] 
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(3.7) 

II  their  is  no  -,uppi  ess  i  on  filter,  hi  )  =  1  for  f  =  ()  and  zero  otherwise. 


There lore,  the  corresponding  output  SNR  is 
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where  o.(0)  represents  the  total  power  of  the  narrow  hand  interference. 

The  ratio  of  the  SNR  in  Hq.  (  3 . i  to  the  SNR  in  liq.  (3.8)  repre¬ 
sents  the  improvement  in  performance  due  to  the  use  of  the  interference 
suppression  filter.  This  ratio,  denoted  by  q,  is 
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We  observe  that  n  is  independent  of  the  processing  gain  1.. 

In  plot  tin  the  improvement  factor,  it  is  convenient  to  use  a 
logarithmic  scale.  Thus,  we  define 

",  lp  =  10  log  n  (3.10) 


This  factor  will  be  plotted  against  the  normalized  SNR  at  the  output  of 
the  I’.N  correlator  when  there  is  no  suppression  filter.  In  other  words, 
the  abscissa  is 


SNR. 


no 


n.(0) 


(3.11) 


As  a  consequence,  the  graphs  of  •  ^  vs.  SNR^/I.  are  universal  plots  in 
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dm 


the  sense  th:it  they  apply  to  any  l’\  spread  spectrum  system  with  arbitran 
processing  gain. 


3.2  Characteristics  of  the  Interference  Suppression  Filter 

In  this  section  we  shall  discuss  the  characteristics  of  the  inter¬ 
ference  suppression  filter  anil  we  shall  illustrate  its  performance  as 
measured  in  terms  of  t he  improvement  factor  n^ 

First  we  consider  two  models  for  the  narrow  band  interference.  One 
type  consists  of  a  sum  of  equally  spaced  sinusoids  covering  20‘j  of  the 
signal  hand. 

In  particular,  i(k)  is  expressed  as 


ilk) 


y 

}  A  cos i It  f 
L ,  m  •  m 
m= 


k  + 


f  ) 

in 


(3.12) 


where  the  amplitudes  (At  were  selected  to  be  identical  and  the  phases 
are  uniformly  distributed  on  the  interval  f  t> ,  2tt  )  .  The  autocorrelation 
I'unct  ion  of  i  (k)  is 


fdk)  =  \  t  A,;  COS  2-  fin  k  (3.13) 

m=  1 

The  number  of  tones  used  in  ilk)  was  either  o  =  100  or  Q  =  10.  A 
second  type  of  narrow  band  interference  employed  was  filtered  white 
noise.  The  filter  characteristic  was  that  of  a  Butte rworth  filter 
having  a  3dB  bandwidth  covering  20'1.  of  the  signal  band. 

figure  3.1  illustrates  the  improvement  factor  for  the  Q  =  100 
sinusoidal  interference  and  the  filtered  white  noise  interference  for  a 
six-pole  Rutterworth  and  a  two-pole  Butterworth  filter.  The  interference 
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Figure  3.1  Improvement  Factor  for  Single-Band  Interference 
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suppression  filter  consists  of  four  taps  with  exact  values  for  the  pre¬ 
dictor  coefficients  determined  from  solving  the  linear  equations  given 
by  (2.17).  We  observe  that  the  improvement  factors  for  the  sinusoidal 
interference  and  the  interference  from  the  six-pole  Buttcrworth  are 
practically  identical.  On  the  other  hand,  the  interference  suppression 
filter  is  not  quite  as  effective  in  suppressing  the  interference  from 
the  two-pole  Butterworth  primarily  because  of  the  difficulty  in  estimating 
and  suppressing  the  rather  significant  amount  of  out-of-band  power. 

Figure  3.2  illustrates  the  improvement  factor  for  the  same  conditions  as 
those  in  Figure  3.1  except  that  the  predictor  coefficients  were  determined 
from  simulation  data  using  the  Levinson  algorithm.  The  results  of  the 
simulation  agreed  very  well  with  the  analytical  results  shown  in  Figure 
5.1  except  at  low  values  of  interference  where  the  improvement  factor 
approaches  zero  dB,  theoretically,  but  the  simulation  data  indicates  a 
small  loss  in  performance.  We  have  observed  this  phenomenon  in  other 
simulation  data,  which  suggests  that  for  small  interference  it  is  best 
to  arbitrarily  set  the  predictor  coefficients  to  zero.  Figure  3.5 
illustrates  the  improvement  factor  for  sinusoidal  interference  with  Q  = 

10  and  Q  =  100  tones.  The  suppression  filter  consisted  of  cither  four 
taps  or  fifteen  taps  with  the  predictor  coefficients  estimated  from 
simulated  data.  There  appears  to  be  little  difference  in  performance 
between  Q  =  10  and  Q  =  100.  Furthermore,  there  is  very  little  gain  in 
performance  when  the  number  of  taps  is  increased  from  four  to  fifteen. 

The  major  conclusion  that  we  have  reached  from  the  above  results  is  that 
the  model  for  the  interference  is  not  critical.  Consequently,  in  most 
of  our  numerical  results  we  used  sinusoidal  interference  with  Q  =  100 
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rigure  3.2  Improvement  Factor  for  Single-Band  Interference 
Simulation  Data 
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Figure  3.3  Improvement  Factor  for  Sinusoidal  Interference 


tones.  To  a  relatively  small  interference  suppression  filter  such 


interference  is  indistinguishable  from  filtered  white  noise. 

Secondly,  we  have  investigated  the  length  of  the  interference  sup¬ 
pression  filter  required  to  achieve  good  performance.  In  this  computation 
we  maintained  the  20“  bandwidth  occupancy  for  the  interference,  hut  we 
distributed  it  equally  in  several  non-overlapping  frequency  hands, 
figure  3.4  illustrates  the  improvement  factor  as  a  function  of  the  number 
of  filter  taps  when  the  SNR  per  chip  without  filtering  is  -20dB.  From 
the  graphs  we  observe  that  a  filter  having  about  eight  taps  performs  well 
when  the  interference  is  split  into  two  frequency  hands;  whereas  a  filter 
having  sixteen  to  eighteen  taps  is  required  to  achieve  good  performance 
when  the  interference  is  split  into  four  frequency  hands.  Figure  3.5 
illustrates  the  performance  as  a  function  of  the  SNR  per  chip  for  eight- 
tap  and  sixteen-tap  filters  when  the  interference  occupies  two  bands  and 
four  hands,  respectively.  The  graphs  show  that  the  sixteen-tap  filter 
with  the  four  interference  hands  closely  approaches  the  performance  gain 
of  the  eight-tap  filter  that  suppresses  the  interference  in  two  hands. 

The  frequency  response  characteristics  of  the  eight -tap  and  sixteen- 
tap  filters  for  an  SNR  per  chip  of  - 20dB  are  shown  in  Figures  5. fa  and  5.7, 
respectively.  It  appears  that  the  sixteen-tap  filter  introduces  some 
distortion  in  the  frequency  range  between  notches.  On  the  other  hand, 
if  the  number  of  taps  is  increased  beyond  sixteen,  the  frequency  response 
is  improved.  For  example,  Figure  3.«S  illustrates  the  frequency  response 
characteristic  of  a  20- tap  filter  when  the  SNR  per  chip  without  filtering 
is  -  2'klB. 

The  computations  reported  above  were  repeated  for  eight-band  and 
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Figure  3.7  frequency  Response  for  Sixteen-Tap  Filter 
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sixteen-hand  interference.  Figure  3.9  shows  the  improvement  factor  as 
a  function  of  the  number  of  filter  taps  when  the  SNR  per  chip  is  -20dB. 
There  appear  to  be  two  threshold  regions  in  these  graphs.  IVhen  the 
interference  occupies  eight  bands,  the  first  threshold  occurs  at  lb 
taps,  where  basically  one  comp  lex -conjugate  zero  pair  is  assigned  to 
each  interference  band.  The  second  threshold  occurs  at  32  taps,  in  which 
case  there  are  two  comp lex- conjugate  zero  pairs  assigned  to  each  inter¬ 
ference  band.  When  the  interference  occupies  sixteen  bands,  the  first 
threshold  occurs  at  .32  taps  and  the  second  threshold  occurs  at  about  b4 
taps.  The  improvement  factor  as  a  function  of  the  SNR  per  chip  without 
filtering  is  shown  in  Figure  3.10  for  a  40-tap  predictor  operating  with 
eight  interference  bands  and  an  80-tap  predictor  operating  with  sixteen 
interference  bands.  The  performance  gain  is  very  similar  to  that  ob¬ 
tained  when  the  interference  is  spread  over  fewer  hands. 

The  conclusion  that  we  have  reached  from  observation  of  the  above 
results  is  that  the  filter  will  suppress  the  multi-band  interference 
provided  that  it  has  enough  degrees  of  freedom,  i.c.,  it  is  sufficiently 
long,  to  assign  at  least  one  complex-conjugate  pair  of  zeros  to  each  band. 
This  behavior  is  substantiated  further  by  the  following  frequency  response 
characteristics  of  the  suppression  filter  when  there  arc  eight  bands  of 
interference.  Figures  .3.11  through  3.14  illustrate  the  frequency  response 
of  filters  corresponding  to  predictors  of  order  8,  lb,  32  and  48,  re¬ 
spectively.  We  observe  that  a  filter  with  eight  degrees  of  freedom  is 
basically  an  all-pass  filter.  It  docs  not  have  enough  degrees  of  freedom 
to  place  notches  at  the  eight  frequency  bands.  On  the  other  hand,  with 
lb  degrees  of  freedom  it  does  place  the  notches  at  the  desired  frequencies. 
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Figure  3 . 10 
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lor  higher- order  predictors,  the  frequency  response  is  improved  as 
shown  in  figures  5.13  and  3.14. 


Another  view  of  the  characteristics  of  the  suppression  filter  is 
provided  by  the  position  of  its  zeros  in  the  --plane.  For  example, 
figure  5.15  illustrates  one-half  of  the  unit  circle  with  the  positions 
of  the  zeros  for  the  suppression  filter  corresponding  to  2,  4,  6,  8, 

10.  12-order  predictors.  In  this  case,  the  interference  is  concentrated 
in  two  hands.  As  shown  in  the  plot,  a  second-order  predictor  places  its 
comp  lex -conjugate  pair  of  zeros  far  from  the  unit  circle  and  roughly 
midway  (in  angle!  between  the  two  interference  hands.  Thus,  its  per- 
tomance  is  poor.  However,  a  fourth  or  higher-order  predictor  does  have 
zeros  within  the  interference  regions. 

The  above  results  indicate  that  a  prediction  filter  having  a  number 
of  coefficients  that  is  fewer  than  twice  the  number  of  interference  bands 
is  useless,  in  the  sense  that  it  does  nothing,  i.e.,  it  is  an  all-pass 
filter.  Apparently,  this  is  a  limitation  of  the  mean  square  error 
criterion  used  to  design  the  prediction  filter.  If  one  knows  that  the 
number  of  degrees  of  freedom  is  fewer  than  twice  the  number  of  inter¬ 
ference  bands,  an  ad  hoc  scheme  such  as  arbitrarily  assigning  a  complex- 
conjugate  pair  of  zeros  to  each  band,  up  to  the  maximum  number  of  bands 
that  can  be  suppressed  with  the  given  number  of  degrees  of  freedom, 
appears  to  he  better. 

Since  we  are  dealing  with  a  digital  communication  problem,  we  have 
also  investigated  the  characteristics  of  the  interference  suppression 
filter  as  a  noise-whitening  filter  in  a  matched  filter  realization.  That 
is,  we  view  the  combined  narrow  hand  interference  plus  wideband  noise  as 
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for  2,  4,  6,  8,  10,  and  12-0 rder  Predictors 


an  cquiva lent  colon.il  noise  process.  Now,  in  detection  ol  a  signal  Slf) 
in  colored  noise  with  power  spectral  density  Pi  i  )  ,  the-  output  SNR  at  the 
receiver  is  maximized  when  the  receiver  consists  of  a  noise-whitening 
filter,  say  11(f),  followed  by  a  filter  matched  to  11(f)  S(f).  Thus,  a 
matched  filter  with  frequency  response  characteristic  !!*(!')  S*(f)  will 
maximize  the  output  SNR,  If  11(f)  represents  the  interference  suppression 
filter  with  impulse  response  hit),  then  ll*( f)  represents  a  filter  with 
impulse  response  hi  - 1 ) .  Thus,  the  cascade  of  these  two  filters  is  a 
filter  having  an  even  impulse  response.  Since  we  have  determined  the 
coefficients  of  11(f)  by  means  of  linear  prediction,  the  coefficients  of 
l!*(f)  are  simply  the  time  reverse  of  those  obtained  for  11(f).  Therefore, 
the  cascade  of  11(f)  and  H*(l)  results  in  a  linear  phase  filter.  Use  of 
such  a  filter  prior  to  the  PN  correlator  improves  performance.  This  is 
illustrated  in  figure  a. lb  for  a  four- tap  and  a  fifteen-tap  predictor. 
iVe  observe  that  at  -20dB  per  chip  SNR  the  four-tap  predictor  in  cascade 
with  its  matched  filter  provides  about  21dB  of  improvement.  The  fifteen- 
tap  predictor  with  its  matched  filter  provides  about  2  IdB  of  improve¬ 
ment.  In  comparison,  the  four-tap  predictor  without  its  matched  filter 
provides  about  l.xlB  of  improvement.  Therefore,  the  inclusion  of  the 
matched  filter  has  resulted  in  about  8dB  gain  at  a  -2(>dB  SNR  per  chip. 

Such  a  large  gain  is  highly  significant  and  suggests  that  the  use  of  the 
matched  filter  is  very  desirable. 

finally,  wv  turn  our  attention  to  the  performance  characteristics 
of  the  filters  designed  from  the  data  by  means  of  the  prediction  algorithms 
described  in  Section  II.  for  this  discussion  it  suffices  to  consider  onl> 
single  band  interference.  The  main  point  that  we  wish  to  make  with  regard 


Figure  3.16  Improvement  Factor  for  Predictor  in  Cascade  with 
its  Matched  Filter 
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to  the  performance  of  tlie  algorithms  is  best  illustrated  with  the  results 
shown  in  figure  3.1'.  We  used  a  block  oi'  fifty  data  samples  to  compute 
the  coefficients  of  a  fourth-order  predictor  by  means  of  the  least 
squares  algorithm,  the  Burg  ajgoritiun  and  the  Levinson  algorithm.  Lor 
the  latter,  the  data  was  used  to  generate  the  estimate  of  the  autocor¬ 
relation  function.  The  predictor  coefficients  obtained  from  the  data 
were  used  in  the  computation  of  the  improvement  factor  given  by  Lq .  (3.9). 
The  graphs  indicate  that  ail  three  algorithms  perform  equally  well.  In 
other  words,  the  difference  in  performance  among  the  three  a  Igor  it  Inns  is 
insigni f icant .  This  behavior  is  furthci  substantiated  by  observing  the 

corresponding  frequency  response  characterist ics  of  the  suppression 

/ 

filter,  lor  example,  figures  3. IS,  .3.19  ;/nd  3. 2d  illustrate  the  fre¬ 
quency  response  characteristics  of  the  suppression  filter  designed  from 
fifty  samples  of  data  on  the  basis  of  the  three  algorithms.  Here,  we 
also  observe  very  minor  differences  in  the  frequency  response  character¬ 
istics.  On  the  other  hand,  when  the  order  of  the  predictor  is  a  large 
fraction  of  the  data  record  length  N,  the  Burg  algorithm  and  the  least 
squares  algorithm  ate  expected  to  yield  better  performance  relative  to 
the  Levinson  algorithm. 

3.3  Perfomvince  of  Interference  Suppress i- m  filter  Based  on  Nonparamctric 

Spect  nil  lyst  i mates 

In  Section.  2.1  uc  described  a  method  for  designing  an  interference 
suppress  ion  filter  based  on  conventional,  nonparametric  methods  for 
spectral  estimation.  ,As  an  illustration  of  the  effectiveness  of  this 
approach,  we  computed  estimates  of  the  power  spectral  density  from  simu¬ 
lated  received  data  and  used  the  resultant  estimates  to  specif)  a  filter 
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characteristic  in  accordance  with  Fq.  (2.2j.  The  helch  method  j  3 1  was 
used  to  generate  the  estimates  of  the  power  spectral  density,  lor  this 
computation,  the  ITT  size  selected  was  6-1  points.  The  resultant  inter¬ 
ference  suppression  filter  consists  of  fifteen  taps.  The  number  of 
data  points  used  to  generate  the  spectral  estimate  is  992. 

Figure  3.21  illustrates  the  estimate  of  the  power  spectral  density 
for  an  SNR  per  chip  of  -20dB.  The  corresponding  frequency  response  of 
the  fifteen-tap  interference  suppression  filter  is  shown  in  figure  3.22. 
The  interference  occupied  20*  of  the  signal  band  as  shown  in  the  graph  of 
the  estimate,  and  the  filter  contains  a  notch  in  the  desired  frequency 
band,  for  comparison,  the  spectral  estimate  shown  in  figure  3.23  is  for 
an  SNR  of  -10dB  per  chip  without  filtering.  The  notch  in  the  filter  is 
now  more  shallow  as  shown  in  figure  3.24.  This  behavior  is  similar  to 
that  obtained  by  means  of  linear  prediction. 

As  a  final  computation,  the  coefficients  h(n)  of  the  interference 
suppression  filter  were  substituted  into  f.q.  (3.9)  and  the  improvement 
factor  was  evaluated,  figure  3.23  illustrates  the  improvement  factor 
as  a  function  of  the  SNR  per  chip  without  filtering.  When  compared  with 
our  previous  results  for  single  band  interference  using  linear  prediction, 
we  find  that  the  performance  improvement  has  similar  characteristics. 
Therefore,  it  appears  that  the  nonpa ramet r i c  method  for  spectral  esti¬ 
mation  coupled  with  the  filter  design  formula  in  liq.  (2.2)  provides  a 
viable  means  for  suppressing  narrow  band  interference  in  a  wideband 
signal.  The  one  disadvantage  of  this  method  is  the  relatively  large 
sample'  sire  required  to  generate  the  spectral  estimate. 
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□  IS  TAP  LIN  PHASE  FILTER 
BASED  ON  WELCH  METHOD 


Figure  3.25  Improvement  Factor  for  15-Tap  Filter 
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5.  I  Per  fonnance  of  Equalized  P.\  Spread  Spectrum  System 

This  section  of  the  report  is  devoted  to  the  per l’o nuance  of  the 
receiver  in  the  presence  of  fading  multipath  signal  components  and  narrow 
band  interference.  As  indicated  previously  in  Section  2.5,  a  major  con¬ 
cern  with  the  use  of  interference  suppression  in  the  presence  of  multipath 
is  the  sensitivit)  of  the  predictor  to  the  channel  multipath  structure. 
Ideally,  one  would  like  to  have  the  prediction  filter  respond  only  to 
the  interference  components.  However,  this  is  possible  only  if  (11  the 
delay  M,  as  illustrated  in  Figure  2.2,  exceeds  the  duration  of  the  channel 
multipath  spread,  i.e.,  use  an  M-step  predictor,  or  (2)  the  delay  M  plus 
the  time  span  of  the  prediction  filter  is  less  than  or  equal  to  the  smal¬ 
lest  time  interval  between  successive  multipath  components  of  the  received 
signal.  The  problem  with  condition  (2)  is  that  it  places  an  unreasonably 
severe  constraint  on  the  length  of  the  prediction  filter.  In  practice, 
the  multipath  characteristics  of  the  channel  are  not  controllable  to  the 
extent  necessary  to  satisfy  condition  (2).  Consequently,  this  condition 
cannot  be  achieved  realistically.  On  the  other  hand,  imposition  of  con¬ 
dition  (1J  results  in  relatively  poor  estimates  of  the  narrow  band  inter¬ 
ference  when  the  multipath  spread  is  large.  This  is  due  to  the  fact  that 
the  interference  i(t)  ;md  its  delayed  version  i(t  -  M)  are  highly  de- 
correlated  when  M  is  large.  Consequently,  the  estimate  of  the  inter¬ 
ference  at  the  output  of  the  prediction  filter  is  poor.  This  phenomenon 
was  observed  in  simulation  results  on  a  two-path  channel  characteristic. 
The  conclusion  is  that  imposition  of  condition  (1)  leads  to  such  poor 
performance  to  render  it  impractical. 

Since  the  M-step  predictor  yields  poor  performance  for  large  M  due 


to  the  decorrelation  in  the  interference  signal,  we  invest  rented  the 
use  of  a  one-step  predictor  in  the  presence  of  multipath.  I  ho  following 
frequency  response  characteristics  of  the  interference  suppression  !  liter 
resulted  from  one- step  prediction  of  a  signal  cons  i.-.  ting  of  narrow  rand 
interference  plus  two  multipath  signal  components  of  equal  strength. 
Single-hand  interference  was  employed  in  this  computation.  The  exact 
autocorrelation  function  was  used  in  the  computation  of  the  coefficients 
of  a  fourth-order  predictor  according  to  l'q.  (d.h). 

figures  3. do  through  5.50  illustrate  the  frequency  response  of  the 
suppression  filter  for  varying  amounts  of  narrow  band  interference  and  a 
multipath  spread  of  two  chips,  lie  observe  that  a  strong  narrow  hand 
interference  completely  masks  the  signal  multipath  components  and,  hence, 
the  filter,  or  one-step  predictor,  responds  to  the  interference  by  placing 
a  notch  at  the  desired  frequency  band.  On  the  other  hand,  when  the  inter 
ference  is  weak,  the  predictor  responds  primarily. to  the  multipath  compo¬ 
nents  and  attempts  to  implement  an  inverse  channel  filter.  The  transition 
point  between  a  strong  interference  and  a  weak  interference  appears  to  Iv 
in  the  range  where  the  signal  and  interference  are  comparable  in  power. 

A  similar  set  of  curves  was  obtained  for  a  multipath  spread  of  four 
chips  as  illustrated  by  the  frequence  responses  in  figures  5.51  through 
5.55.  Here,  again,  the  suppression  filter  adapt.--  to  the'  inverse  channel 
filter  when  the  interference  i~-  small,  but  is  unaffected  f\  the  multipath 
when  the  interfeience  is  large.  On  the  other  hand,  if  we  increase  the 
multipath  spread  beyond  four  chips,  the  four-tap  predictor  will  not  see 
the  multipath  component.  This  condition  is  illustrated  by  the  frequence 
res pon .-e  curve,  m  1  igures  5.5o  through  5.10  for  a  multipath  spread 
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Figure  3.32  Frequency  Response  of  4-Tap  Filter  for  Multipath  Spread 
of  Four  Chips 
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corresponding  to  five  chips.  In  this  case,  when  the  interference  is 
small,  the  suppression  filter  approaches  an  all  pass  filter. 

The  results  given  above  demonstrate  that  the  one-step  predictor 
can  be  used  effectively  in  a  multipath  environment.  On  the  other  hand, 
the  M-  step  predictor,  vdien  M  is  large,  is  not  very  effect  ive  except , 
perhaps,  for  a  pure  vU  interference. 

An  important  issue  in  our  investigation  of  narrow  band  interference 
suppression  in  a  l’.\  spread  spectrum  system  is  the  amount  of  intersymbol 
interference  introduced  by  the  suppression  filter  and  whether  or  not 
there  is  a  need  for  ei|un !  fiat  ion.  In  order  to  demonstrate  tin-  magi: I  rude 
of  the  distortion  »,  uused  bv  lie.  interu  rence  suppression  filter,  we  have 
plotted  in  figure  3.11  rlie  variance  of  each  of  the  three  terms  in  the 
denominator  of  The  expression  tor  tin.  SNR  given  by  fq.  (3.').  In  this 
c-Huput  .it  ion.  the  channel  is  ideal  i  no  multipath)  and  the  interference  is 
confined  to  a  single  fund.  The  suppression  filter  consists  of  the  filter 
A,J~l  in  cascade  with  its  matched  filter.  The  term  "self-noise"  refers 
to  the  time  dispersive  distortion  caused  bv  the  suppression  filter.  It 
is  evident  that  this  is  the  dominant  term  in  the  denominator  of  the  ex¬ 
pression  for  the  S\  R  in  fq.  l3.'-i.  I  lent,  e,  it  is  the  term  that  limits 
the  performance  of  the  svstem  when  the  narrow  hand  interference  is  strong, 
for  example,  at  an  SNR  per  chip  of  -.llhlf.  the  distortion  due  to  the  filter 
i  .  mine  than  I3.IB  above  the  white  noise  and  more  than  Jo. If  above  the 
residual  int  erlcivnce.  Alien  the  interference  is  weak,  the  filter  approache: 
an  all -pass  eharactei ist  ic  and.  hence,  the  distortion  due  to  time  dis 
per-i.n  >  it  c  rt  a-c-  .  Therefore,  the  need  for  equalisation  of  the  time- 
dispersive  et  iVcts  of  the  suppression  filter  is  established  by  these 
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1  he  results  of  the  Monte  Carlo  simulations  are  shown  in  Figures 
3.1.  through  3.4b.  The  first  two  figures  illustrate  the  error  probability 
for  a  signal -to- interference  ratio  of  -20dB  ;tnd  -lOdR  per  chip  prior  to 


filtering,  respectively.  In  this  case,  the  input  to  the  PN  correlator 
consists  of  the  chip  estimates  from  the  PFH  prior  to  binary'  quantitation, 
figure  3.44  illustrates  the  effectiveless  of  the  interference  suppression 
filter.  That  is,  the  performance  of  the  system  with  a  signal -to- inter¬ 
ference  ratio  of  -20dB  is  almost  as  good  as  that  for  -lOdB.  A  check  of 
the  slope  of  these  error  rate  curves  indicates  that  they  correspond  to 
the  performance  of  a  dual  diversity  system.  In  other  words,  the  feed¬ 
forward  part  of  the  equalizer  acts  as  an  equivalent  (coherent)  diversity 
combiner,  figures  3. 43  and  3.4o  illustrate  the  performance  of  the  re¬ 
ceiver  when  the  input  to  the  PN  correlator  are  the  hard  decisions  from 
the  DM:.  Comparison  of  these  graphs  with  those  in  Figures  3.42  and  3.43 
leads  us  to  conclude  that  chip  estimates  into  the  PN  correlator  result  in 
better  performance  at  low  error  rates. 

Fhe  error  rate  results  given  above  indicate  that  inter-chip  inter¬ 
ference  caused  by  channel  multipath  and  the  suppression  filter  can  be 
compensated  by  the  DIF.  It  is  important  to  emphasize,  however,  that,  for 
the  Dll:  to  be  effective,  it  must  be  sufficiently  long  to  span  the  time- 
dispersion  of  the  filter  anti  the  channel  multipath.  That  is,  an  increase 
in  the  length  of  the  interference  suppression  filter  must  be  accompanied 
by  an  increase  in  the  length  of  the  Dll'. 
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MISSION 

of 

Rome  Air  Development  Center 

RA VC  plan*  and  execute*  research,  development,  test  and 
selected  acquisition  programs  In  support  of  Command,  Control 
Communications  and  Intelligence  (C3I)  activities .  Technical 
and  engineering  support  utithln  areas  of  technical  competence 
Is  provided  to  ESP  Program  Offices  ( P0-& )  and  other  ESP 
elements.  The  principal  technical  mission  areas  are 
communications,  electromagnetic  guidance  and  control,  sur¬ 
veillance  of  ground  and  aerospace  objects.  Intelligence  data 
collection  and  handling,  Information  system  technology. 
Ionospheric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  maintainability  and 
compatibility. 
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